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Pion in the Holographic Model with 5D Yang-Mills Fields 
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We study pion in the holographic model of Hirn and Sanz which contains two Yang-Mills fields 
defined in the background of the sliced AdS space. The infrared boundary conditions imposed on 
these fields generate the spontaneous breaking of the chiral symmetry down to its vector subgroup. 
Within the framework of this model, we get an analytic expression for the pion form factor and a 
compact result for its radius. We also extend the holographic model to include Chern-Simons term 
which is required to reproduce the appropriate axial anomaly of QCD. As a result, we calculate 
the anomalous form factor of the pion and predict its Q 2 -slope for the kinematics when one of the 
photons is almost on-shell. We also observe that the anomalous form factor with one real and one 
virtual photon is given by the same analytic expression as the electromagnetic form factor of a 
charged pion. One of the advantages of the present model is that it does not require an infrared 
boundary counterterm to correctly reproduce the anomaly of QCD. 

PACS numbers: 11.25.Tq, ll.10.Kk, 11.15.Tk 12.38.Lg 



I. INTRODUCTION 

During the last few years applications of gauge/gravity 
duality [1] to hadronic physics attracted a lot of atten- 
tion, and various holographic dual models of QCD were 
proposed in the literature (see, e.g., [l]-[23j]). These 
models were able to incorporate such essential proper- 
ties of QCD as confinement and dynamical chiral sym- 
metry breaking, and also to reproduce many of the static 
hadronic observables (decay constants, masses), with val- 
ues rather close to the experimental ones. 

In our recent papers [H, [2{| [2(| we developed a formal- 
ism that allows to systematically study the meson form 
factors within the holographic "hard-wall" approach of 
Refs. @,[|. We applied it first to form factors and wave 
functions of vector mesons 0, [H| and the n l26ll to the 
pion electromagnetic form factor. In Ref. [231, we ex- 
tended the holographic dual model of QCD to incorpo- 
rate the anomalous F 1 , 1 *^o(Q\,Q\) form factor. 

In the present paper, we consider a holographic model 
of QCD proposed by Hirn and Sanz pjj], with SU(2)l x 
SU{2) R Yang-Mills (YM) gauge fields living in the back- 
ground of sliced five-dimensional (5D) AdS space. Un- 
like the approach of Refs. [f| Q, this model does not 
require the existence of an additional degree of freedom 
dual to the chiral condensate of four-dimensional (4D) 
QCD, which spontaneously breaks the chiral symmetry 
via the Higgs-like mechanism. Instead, the chiral symme- 
try breaking down to its vector subgroup SU(2)y occurs 
due to the boundary conditions (b.c.) imposed on the 
infrared (IR) brane. 

At the same time, the global SU (2) r, x SU (2)r symme- 
try of QCD is generated from the requirement that the 
fields vanish on the ultraviolet (UV) boundary. The chi- 
ral field U(x), the phase of which describes the Nambu- 
Goldstone bosons, appears in this model as a product of 
Wilson lines connecting IR and UV branes by the fifth 
components of left and right gauge fields. 



Since the model of Ref. [l2| incorporates YM fields 
only, it has a single parameter, the cutoff scale zq, that 
determines the size of masses m p , m ai , etc., of vector and 
axial- vector mesons, and also the pion decay constant / ff . 
Still, many features of this construction are similar to 
those discussed in the papers @, @, [28| . 

The paper is organized in the following way. We start 
by outlining, in Section II, the basics of the hard-wall 
model of Hirn and Sanz. In particular, we write the form 
of the 5D action, show how to separate gauge fields into 
dynamical and source parts, define the chiral fields as 
Wilson lines, and demonstrate how the boundary condi- 
tions on the fields break the global symmetry of QCD 
down to the vector subgroup. We also elaborate on the 
meaning of the boundary conditions, and present vector 
and axial-vector fields in a way helpful for further studies. 

In Section III, we calculate and analyze the pion form 
factor. We show that the form factor can be represented 
analytically in terms of the modified Bessel function, and 
obtain a compact analytic result for the pion charge ra- 
dius. We further explore the large-Q 2 behavior of the 
form factor and observe good agreement with experimen- 
tal data. 

In Section IV, we consider the generalization of 
the AdS/QCD model that includes isoscalar fields and 
Chern-Simons term. Using this extended model we de- 
scribe the calculation of the F 7 . 7 » 7r o(Q 2 , Q|) form factor 
and express it in terms of the pion wave function and two 
bulk-to-boundary propagators for the vector currents de- 
scribing EM sources. We observe that in case of one real 
photon, the anomalous form factor of the neutral pion 
is identical to the electromagnetic form factor of charged 
pion. We discuss kinematics with one real and one vir- 
tual photon and calculate the value of the <5 2 -slope of the 
form factor. We also investigate the formal limit of large 
photon virtualities, and compare these results to those 
obtained in our earlier paper [26( . Finally, we summarize 
the paper. 
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II. OUTLINE OF THE MODEL 
A. The Setup 

The model of Ref. 1131 is based on the action 



Sym = ~T~2 J d 5 x^/g Tr 



L 



AIN 



RmnR 



MN 



(1) 



with the metric 



g M Ndx M dx N = — (rj^dx^dx" - dz 2 ) , (2) 

where z e (0, z ], r)^ = Diag(l, — 1, — 1, — 1), (i, v = 
(0,1,2,3), M,N = (0,1,2,3,*), 

A M n = d M A N - d N A M - i[A M ,A N ] , (3) 

and A M = t a A a M , A = {L, R}, (t a = a a /2, with a a being 
Pauli matrices). The gauge fields transform as 

A M {x, z) -> g A A M g^(x, z) + igAdM9~^(x, z) , (4) 

where gA (x, z) £ SU(2)a- On the UV brane, the bound- 
ary conditions L^x, 0) = t^x) and 0) = r^[x) are 
assumed, where ^(x) and r^(x) are the sources for the 
left and right 4D currents. Vector = (L^ + i? M )/2 
and axial- vector — (L^ — R^)/2 gauge fields are 
dual to the vector and axial-vector currents of QCD re- 
spectively. Working in the axial-like gauge, in which 
L z (x, z) = R z (x,z) = 0, one can write the vector 
and the axial-vector A^ fields as 

% (x, z) = (x, z) + (x, 0) , (5) 
A^ (x, z) = A^ [x, z) + a (z) A^ (x, 0) , 

where the so called "dynamical" fields (x, z) and 
A^ [x, z) satisfy the following b.c. 



U M (x,0) = 0, A^x,0) = Q 



(6) 



on the UV brane. However, on the IR brane, the vector 
field obeys Neumann b.c. 



d z V,j,{x,zo) = , 



(7) 



while both of the axial-vector fields A^ and A^ are re- 



quired to satisfy Dirichlet b.c. 



Ap(x, z ) = , A^ (x, z ) = 



(8) 



As pointed out in Ref. [l2| (and will be discussed be- 
low), in order to avoid the mixing between the pion and 
the axial resonances, the function a (z) should satisfy the 
equation 



d z (^g^g zz d z a{z)) = . 



(9) 



The following b.c. on the function a(z) 

a(0) = l, a(z o )=0, (10) 
are determined from the b.c. in ^ and ([5]). As a result, 



>,(z) = l-z 2 /zl 



11) 



The chiral field 

U(x)=Z R (x)il 1 {x) i (12) 
is built from the path-ordered Wilson lines: 

£l(x) — Pexp i-i J dz'L z (x,z')^ , (13) 

£,r{x) — Fexp |— i J dz'R z (x, z') | . 

With respect to the global chiral transformations, the 
field U(x) transforms in the same way as the chiral field 
in the non-linear sigma model. Therefore, the pion field 
is build from a product of Wilson lines extending from 
one boundary to the other. 

B. Meaning of Boundary Conditions 

The Dirichlet b.c. V^O) = and A^(x,0) = 0, 
imposed on the UV brane, are equivalent to L^(a;,0) — 
R^(x,0) = 0. The latter b.c. are important, since for 
these the residual gauge invariance is a global symmetry 
of 4D QCD. Another significance of these b.c. is that 
they secure a finite action at the UV boundary. Indeed, 
since the lagrangian ^/gTrF 2 is singular at z = 0, the 
field strengths have to vanish there to produce a finite 
action. We can partially fix the gauge by requiring that 



L M (x,z^ 0) = , R M (x,z^0) =0 



(14) 



This gauge choice remains unaltered if we perform ad- 
ditional gauge transformations g les (x,z — > 0) satisfying 
the condition 



d M g TCS (x,z^0) = , 



(15) 



which means that g ICS (x,z) goes to a constant matrix 
9l,r S SU(2)l,r at z = 0. In the holographic model, 
{9l, gn) € SU(2)l x SU(2)r corresponds to the global 
chiral symmetry of QCD, at z = 0. 

The other Dirichlet b.c. A fJ _(x, zq) = (imposed on 
the IR brane) breaks gauge invariance in the bulk, re- 
quiring Lfj,(x, zq) — R^(x, zq), which is equivalent to the 
condition g r £ s (x,z — z ) — g r ^ s (x,z = z ) = h(x). The 
resulting breaking of gauge invariance in the bulk leads 
to the spontaneous breaking of the global chiral symme- 
try on the 4D UV brane down to the vector subgroup. 
As a consequence, the Wilson lines £,l,r transform as 



£,l,r = 9l,r £,l,r h[x)^ , 



(16) 
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where h(x) G SU(2)v is a local gauge symmetry on the 
4D IR brane. 

Finally, the remaining Neumann (gauge invariant) b.c. 
V Z n(x, Zq) — d z Vfj,(x, zq) — is required to have a unique 
solution for the equations of motion. This b.c. was cho- 
sen to preserve the vector gauge invariance, since other- 
wise the breaking of it may lead to spontaneous breaking 
of the global vector symmetry of QCD. However, accord- 
ing to the Vafa-Witten theorem [29( this does not occur 
in QCD. 



C. Generalities 

It is useful to define the following 4D fields: 
u fl (x) = i{^ R (d fl -ir fi )^ R -^(d fl -U IM )^ L ] , (17) 

r„ (*) = \ {in (d„ - ir M ) U + & (0„ - U„) £ L } . 
Notice that 



L zli = ^ L ( d z V^ + 8 Z A^ - - (d z a) Uli )e L , (18) 



R 



- £,r \ d x Vn - d z A fi + - (d z a) u M ) £ R , 



where 



Za{x, z) = Pexp < — i / dz'A z (x, z') 



(19) 



It is straightforward to see that Z,a{x,z = 0) = £a{x). 
The field T M (x) will be used later to define the covariant 
derivative. 

Now, we are in a position to explain the particular 
choice of the field a(z) described earlier. Indeed, accord- 
ing to Eq. (fl8|). the term in the 5D action (|4ip . which 
describes the mixing of axial-vector and pion fields, is 
proportional to the integral 



z ° dz 

z 



P z ° dz ( 1 \ 
(d z AJ(d z a) = - J —AJd z -d z a\ . (20) 



In order to avoid this mixing, one may impose the EOM 
© for the field a(z), namely, 



•->-. I - d z a I =1) 



(21) 



As a result, the integral in Eq. (f2TJ)) vanishes automati- 
cally. 

It is instructive to observe that the part of the YM 
action with 4D indices can be written as 



1 



Tr(Rl v + Ll v ) = ~Tr(Fl tl „+F^ v 



where 



F±^ (x, z) = SlL^Zi. ± ZrR^R • (22) 



The field strength tensors of the sources are defined by 

f ±ttv (x) = zU^Zl ± • (23) 

One can rewrite the field F +flL , as follows 

F+tiu = 2 (V M K - V„y M - i [V^ V u ] - i [Afj,, A,]) (24) 

1 - a 2 

+ ia ([u^Au] + [AftjU,,]) + + i — - — [u^Uv] , 

where the covariant derivatives of the fields are given by 
V M - = ■ + [r M , •]. In the same way, the field F_ fll/ can 
be rewritten as 

F^ v = 2 (V M A„ - - i [V^A U ] - i [Ap, V v \) 

+ ia ([«„, V v ] + \y M , u v ]) + af-^ . (25) 

Notice, that the fields A^ and are dynamical fields (to 
be discussed in details in the sections below) . These fields 
contain axial- vector and vector mesons only. Information 
about the pion field is contained in the fields m m and T^. 

D. Vector fields 

The dynamical vector fields have the following repre- 
sentation 



V„(x,z) = J2v<r\x)il> n (z) , 



(26) 



n=l 



in terms of the wave functions ip n (z) satisfying EOM 

[z 2 d 2 - zd z + M 2 z 2 ] ip n {z) = , (27) 

with b.c. Vn(0) = d z ip n (z Q ) — 0. Here, e.g. the field 
V^ 1 \x) = g§Pn{x) describes the p-meson. The solution 

for ip n (z) is 



V2 



ZoJl("yo,n) 



zJi{M n z) 



(28) 



where M n is determined from Jo(M n Zo) = and, there- 
fore, M n = 7o : n/^o (with J (7o : n) = 0). The value of 
zo = 1/ (323 MeV) is fixed from the experimental mass 
of the p-meson Mi = 776 MeV. Eigenfunctions ip n are 
normalized as 



dz 



(29) 



The Fourier transform of the vector field is written 
as Vfjt(q, z) — V/j,(q)V(q, z), where V^{q) is the Fourier 
transform of the 4-dimensional field V^x), and V(q,z) 
is the bulk-to-boundary propagator. The latter satisfies 
the EOM 

zd z (^d z V(q,z)\ +q 2 V(q,z) = (30) 
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with b.c. V(q, 0) = 1 and d z V(q, zq) — 0. It can be also 
written as the sum 



- Ml 



(31) 



where f n is the decay constant of n th vector meson and 



fn 



-d z lp n (z) 



V2M n 



=0 g5ZoJl(lO,n) 



(32) 



E. Axial-vector and pion fields 

The dynamical axial-vector fields can be written as: 

oo 

A fl (x,z) = J2 A( ? ) (x)Tpn(z) , (33) 

n=l 

where the functions Vn( z ) satisfy the same EOM as 
4>n(z), but with different b.c. 4>n (0) — ( z o) — 0- Here, 
in particular, the field A^\x) = g^ai^x) describes a\- 
meson. 

The solution for the axial- vector sector is 

^(z)^z.h(M A z) , 

where is determined from IR b.c. Ji(M^zo) = 0. 

In the axial gauge, the axial-vector field with the dy- 
namical fields turned off is given by 

(x, z) — a (z) A^ (x, 0) 
ia (z 



{&8vt L -&d lt { R } . (34) 



Taking into account the definition of Wilson lines 
^l,r{x), one can check that 



= -(d M tl)tL = ~i \ dz'd^L z {x, z') , 
Jo 



(35) 



and, therefore, 



A a (x, z) = a{z)d li / dz' A a z (x, z') = a(z)(d ll Tr a ) (36) 
Jo 

Notice, that the same result could be obtained even sim- 
pler if one uses the additional gauge redundancy by fixing 
£]r = U and £r = 1, in which case 



(37) 



and since U = e 2l7r , then Ud^U^ = — 2z<9 M 7r, therefore, 
A^ = a(3 M 7r). 



If £ M = r M = and A^ = = 0, one has 



(x) = i{&d^R-&d lt Z I ,} , 



(38) 



L zt _i — 



4 (d z a)£ L u^{ , 



(39) 



Notice, that 



^R U ^R 



(40) 



The order 0(p 2 ) kinetic term in the action for the chiral 
fields U(x) is coming from the following part of the 5D 
action: 



'kin 



1 



d 4 x 



dz 



— Tr (L% + R 2 Z 



(41) 



Taking into account Eqs. (j3"5|) and (|40[) . and integrating 
over z, one obtains that the kinetic term in the action for 
the chiral fields becomes 



S kin = J d 4 x L Tr {d^d^U) 



where 



J 7T 



95 Jo 



dz 2 2 
— (d z a) = -5-5 

z 914 



(42) 



(43) 



Furthermore, integrating by parts and using Eq. (121 
gives 



, 2 1 f*° dz a 1 a{z) 
9%Jo z g*\ z 



z=0 



(44) 



Since a(0) is normalized to 1, the pion decay constant 
/tt is determined by the value of the function a'(z)/z at 
z = 0. This result is similar to that obtained within the 
holographic model of Refs. 0, where f% is given by 
the z — value of the function 



1 1 1 



(45) 



with &(z) being the pion wave function of that model. As 
we argued in Ref. [2y|, it is the function <5(z) ~ ^'{z)/z 
that is the most direct analog of quantum-mechanical 
wave functions of bound states. Thus, in the present 
model we can introduce an analogous function 



<p{z) 



9 r ofir 



■ d z a(z) 



(46) 



which has <p(0) = g^f-n normalization at the origin. 
In fact, given the explicit form of <x[z), one finds that 
ip(z) — g^f-n for all < z < zq. 

Finally, the full axial- vector field in the axial gauge is: 

oo 

A^x,z) = a(z)d^(x) + . (47) 

n=l 
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The longitudinal part of the axial- vector field Ajf (a;, z) 
9 M (7r a a) can be written as 



A\^P,z)=ip^ a {p)a{z) 



(48) 



where A?, (p, z) and 7r a (p) are the Fourier transforms of 
A^(x,z) and Tt a (x), respectively. Furthermore, since 

A \\^ z ) = %(?)«(*)> thcn 



/( P Wz) = -|i;»«(z) . 



This allows us to rewrite Ajj^(p, z) in the form 

4( P ,z)^I; Q ( P ) a (z) 



(49) 



(50) 



involving the longitudinal projector p a p^/p 2 and the pion 
"wave function" a(z). 



F. Two-Point Function 

The spectral representation for the two-point function 
of axial-vector currents can be written as 



2 pV 



t/3 p a p 13 
9 



E 



F A,n 



p 2 -M* A , 



(51) 

(nonpole terms), 



where (0\J%\tt(p)) = if,P a and (0|J^|A„(p, s)) = 
i^ in e"(p, s), -Fa,™ correspond to the n th axial- vector me- 
son decay constant. Finally, as was shown above, the 
pion decay constant in this model is given by 



f 2 = — — 

U n 2 7 2 ' 
95 Z Q 



(52) 



The hard- wall scale is usually fixed from fitting the phys- 
ical mass of the p-meson, which gives zq = 1/(323 MeV). 
The constant g§ is fixed from correspondence between 
AdS/QCD results and the asymptotic behavior of (per- 
turbative) QCD at large Q 2 , in which case g 2 = 6tt 2 /N c , 
and, therefore, 



f = 



3tt 2 z 2 



(53) 



For N c = 3, this gives f n ~ 102.8 MeV instead of /° xp = 
130.7 MeV. Since f 2 ~ 0(N C ), one may speculate that 
the difference between the two values is related to 0(N®) 
corrections to the f 2 prediction of the AdS/QCD model. 



III. PION ELECTROMAGNETIC FORM 
FACTOR 

Our next step is to apply the Hirn-Sanz holographic 
model [l3| , which does not have additional bulk field dual 



to the chiral condensate, for calculation of the pion form 
factor. A brief discussion of this form factor was given in 
the original paper [e| by Hirn and Sanz. Our goal is to 
incorporate the formalism developed in Ref. [26| (see also 
[3l|). where it was applied to the pion form factor within 
the framework of the AdS/QCD model of Refs. @, i]. 



A. Three-point function 



To find the pion form factor, we need to consider three- 
point correlation function between EM current (0) and 

j4^. c ipixi-ip2X2 (54) 



two axial currents J^ a (xi), J§t{x2) 

T IMX p(pi,P2) = J (Pxx J d 4 x 2 e" 

x (0|T^(x 2 )J*(0)J B a(a:i)|0) , 

where pi, p 2 are the corresponding momenta, with the 
momentum transfer carried by the EM source being q = 
P2 — Pi (q 2 = — Q 2 < 0). The spectral representation for 
the three-point function is 

T» af) (px,p2) =pfpg(pi F ff ] (55) 

P1P2 

+ (transverse terms) + (nonpole terms), 

where the pion electromagnetic form factor F T (Q 2 ) is 
defined as 

Mpi)|^(0)Kp 2 )) = F„(q 2 )( Pl + p 2 ) M . (56) 

According to the prescription of Ref. [26j , the pion form 
factor can be obtained from the three-point function us- 
ing 

Pl a p 2 /3T^V,p 2 )| p?=0ip|=0 = (px + P 2rf 2 F„(Q 2 ) . 

(57) 

The part of the 5D lagrangian, which may contribute 
to the pion form factor, is given by 

V?£vm = t4z Tr (R^R^ + L^L^) 



±9p 
1 

W z 



D 



Tr (L% + R 2 zfi ) 



(l-a 2 )Tr(V^ [u^u v ]) 



1 2 

— (d z a) Tru M u M , 



where V^ v = d^V" - d u V^. Taking into account that 
A^y(a;, 0) C — u°/2 (ignoring the sources and r M ), one 
can derive that 



1 



-^(l-a*)e a » c {d»V»> a )A^ A C H 
7 V r n;(9,a) J (i t %) 



6 



The second term was left "as is" , since the sources in , 
in combination with Wilson lines, will also contribute to 
the pion form factor. 

To calculate the 3-point function, we perform first the 
Fourier transformation, so that V"(q,z) — V£(q)V(q,z) 
is an image of V^(x,z), and A°ii(p, 0) is an image of 
.A^ii (x, 0). Then, varying the action corresponding to the 
first term in £ p7r7T with respect to the sources, V"(q), 
A b a ^(pi) and A^A—pi), produces the following 3-point 
function: 

WJv)J\\aM) J uJ-P^ = i(27r)V 4) (<z +pi -P2) 

v a n 13 1 f z " 1 

X e a bc (Pi + PlYiTi q 2 dz- V(q, z) (1 - a 2 ) , 



P1P2 



2.9. 



5 Je 



where, anticipating the limit p\ — > 0, p\ — > 0, we took 
(Pi<z) = —(P2q) = —q 2 /2 in the numerator. We also took 
into account that 



^l|(P»°) = ^^l|(P) 



(58) 



where D£/ = 9C/ + iUi — irU . Expanding U in powers of 
7r produces the local term 



f 2 



Tr (DHJDptf) D fle abc V^{x)TT b {x)d^ c {x) 



that compensates "—1" in Eq. (j6"2"|) . Namely, performing 
Fourier transform, taking into account that 

7r -(p)=-^»| | (p) , 

and V^a:) = [^(^) + r ( a; )]/2, with further varying the 
corresponding action, gives the following result for the 
total pion form factor: 



F^Q 2 



1 



9^ f 2 Jo 



dz z J{Q, z) 



dz zJ(Q, z) <p 2 (z) . 



(63) 



Using explicit form of (p(z) and incorporating the result 
g%fl = 2/z 2 (12) for /„ we obtain 



B. Form factor derivation 

Now, representing 

= i{2v f5^ (q+ Pl -p 2 ) e abc T^ 3 ( Pl , pa) 

and applying the projection suggested by Eq. ([57]), we 
have 

1 r Zt > A-/ 

f > F U(Q 2 ) = -^Q 2 -J(Q,z)[l-a*(z)] , 
z Ss Jo z 

(60) 

where J(Q, z) = V(iQ, z) is the bulk-to-boundary prop- 
agator taken for spacelike momenta and explicitly given 
by 



J(Q,z) = Qz 



K 1 (Qz) + h(Qz 



K (Qz ) 



Io{Qzo) 



(61) 



Integrating by parts Eq. (|60p and using equations of mo- 
tion both for J and a produces 



FA\Q 2 ) 



1 



r *° f d z OL 

2 , 2 / dzzJ(Q,z) 



9b f 2 Jo 



z 



(62) 



FAQ 2 ) = — / dzzJ{Q,z) . (64) 

z JO 

Notice that, since J(0 : z) = 1, we have correct normal- 
ization for the pion form factor F v (0) = 1 . 



C. Results 

Using EOM for J in the last equation and integrating 
further by parts gives 



FAQ 2 ) 



Z /I 

2r ' 1 



, , dzz'd, [ -d z J{Q,z) ) (65) 
Q z o Jo 



1 - 



Q 2 4 



Io(Qzo) 



From the analytic expression for the form factor, it is 
straightforward to obtain the pion electric charge radius: 



(rl) 



-6 



dFAQ 2 ) 
dQ 2 



Qz 2 



Q 2 =0 



Taking zo = 1/(323 Me V) gives numerically 
(r 2 ) ~ 0.42 fm 2 , which may be compared with the 
experimental value 0.45 fm 2 [331 ]. 

In the large-Q 2 limit, it follows from Eq. (|6"5")) that 



Integrating the second term in £ p7r7r with respect to z 
gives: 



Q 2 FAQ 2 ) -> — - 0.42 GeV 2 



(67) 



1 



dz 



4ffg Jo z 



{d z af Tr = f -f Tr [(D"U) (D^)] 



It is interesting to note that the highest Jefferson 
Lab's experimental points correspond to Q 2 F° xp (Q 2 ) ~ 



7 



Q 2 F„(Q 2 ) (GeV 2 ) 
0.5 




FIG. 1: Pion form factor Q 2 F^(Q 2 ) from the holographic 
model (solid curve, red online) in comparison with the 
monopole data fit Q 2 F™ ono {Q 2 ) = Q 2 /(l + Q 2 /0.50GeV 2 ) 
(dashed curve, blue online). 



can extract the pmr coupling: 

<W = T- lim (Q 2 +M 2 )7V(Q 2 ) (69) 

= T y° \ I dxx 2 Ji (70,1a;) . 
Ji(7o,i) Jo 

Here we took into account the correct normalization of 
currents, giving a factor of \[2. From this result it follows 
that numerically = 1.383 35. Taking g 5 = V% 7r we 
obtain g p7r7r ~ 6.15. The experimental value is g p7r -rr — 
6.03 ± 0.07. The numbers obtained in Models A and B 
of Ref. [6j are 4.48 and 5.29 respectively. 



IV. ANOMALY 
A. Chern-Simons action 



0.4 GeV 2 , which is very close to the holographic model 
result of Eq. (|6"7) (see also Fig.l). 

It is also instructive to use Eq. (|43|) to substitute l/z 2 
in terms of f 2 . This gives 



Q 2 F^Q 2 



(68) 



the expression analytically coinciding with our result ob- 
tained in Ref. [26] within the framework of the AdS/QCD 
model of Refs.[i,[|. 

This outcome has very basic reasons. Namely, the 
large-Q 2 behavior of the form factor is determined, first, 
by the large-Q 2 form of the bulk-to-boundary propaga- 
tor J(Q, z) (which coincides with its free-field version 
K,(Qz) = QzK\{Qz) in any model) and, second, by 
the value of the pion wave function ip(z) (or $(2;)) at 
the origin. The latter also determines / ff in both holo- 
graphic models, which results in the same analytic re- 
sult 2g|/ 2 for the large-Q 2 limit of Q 2 F 7r (Q 2 ) when it is 
expressed in terms of /„.. However, as we already dis- 
cussed in Ref. 26], if we take the experimental value 
/° xp ~ 131 MeV, Eq. (JB5J gives the value 0.67 GeV 2 that 
is well above Jefferson Lab's experimental points. In this 
sense, the expression (IBTl) for the limit of Q 2 F 7T (Q 2 
terms of zq is numerically more successful than Eq, 
One may speculate that since the pion form factor in our 
calculation is given by the ratio of the 3-point function 
term (which is proportional to f 2 F n (Q 2 )) to the 2-point 
function term (which is proportional to / 2 ), the overall 
error of the model in the value of f v is cancelled, and the 
remaining expression for F^^Q 2 ) in terms of zq correctly 
reflects information about the pion size. 



From the form factor expression (|63p and the decom- 
position of c7(Q, z) over the Q-channel bound states, one 



Since the Chern-Simons/Wess-Zumino-Witten term 
for 5/7(2) gauge/global group is vanishing, we need to 
extend the flavor symmetry to U(2)l x U(2)r, so that 
the fields are written as 



f a B a 



(70) 



In order not to confuse the hats on the £7(1) fields with 
the hats on the gauge fields in the axial gauge, we will 
assign hats only to the former and omit these for the 
latter. 

Extending the holographic dictionary, the 4D isosinglet 
vector jj 1 I ~ ^(x) current will correspond to 



J^ 1 0} = 2 ( u 7m w + d ln d ) = t;97m 1 <7 -> %{ x > z ) 



(71) 



where V^x, z) e U(l)v is the abelian part of the U{2) v 
field. We also remind that the third component of the 
isovector J^ /_1 ^' a (x) current corresponds to 



Ji 1=lh3 = I 



(u-y^u - dry^d) = qi^—q 



in Note that the EM current of QCD is defined as 



J. 



EM 



— T{I=n,3 + I j{I=0} 
~ J M 3 ^ 



(72) 



(73) 



It has both isovector ( "p-type" ) and isosinglet ( "cj-type" ) 
terms. 

The 0(B 3 ) part of the 5D CS action, in the axial gauge 
B z = is 



M f 
S ™ [B] = 24^ e " VP ° Tr J dAxdz (W 



l>CT 
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where = d^B u — d^B^. In the holographic model (cf. 
[HI), the CS term is 



S$t[B L ,B R ] = S§[B L ] S%[B R ] , (74) 

where B L ,R = V ± A and V(A) € U(2) V{A) . 

After long, but straightforward calculations, we get 



(3)r 



field dual to the chiral condensate operator of QCD, and, 
moreover, we do not need to add a counterterm on the 
IR boundary to reproduce the correct normalization of 
the anomalous form factor. 

When only one of the photons is virtual, Q\ = Q 2 , 
while another is real, Q\ = 0, we have 



5«AdS = __c_^pa I d 



12tt 2 



fZ 

'.C I dz TT a 

'o 



(75) 



2 (d z a) (d p v*) [d a V v \ - ad z (d p v^ d a v t 



Integrating the second term by parts with respect to z 
and taking appropriate care on the IR boundary gives 



~AdS 
5 CS 



N, 



£_ fivpcr 



4tt 2 



dz ( 



(d z a) I d*x n a (d p V;) (d a %) 



Recall that a(z) = 1 — z 2 / z 2 and, in this model, it has 
the meaning of the pion "wave function" . 



B. Anomalous Form Factor 

In QCD, the 7*7*7r form factor is defined by 



d 4 x e-^>,p|T{j£ M (x) J£ M (0)}|0) (77) 

= e^ af} q la q 2f3 F 1 , 1 ^ (QlQ 2 ) , 

where p = q\ + q 2 is the pion momentum, q± , q 2 are the 
momenta of photons, and q\ 2 = —Q 2 2 - 
Varying Sqq S gives the 3-point function: 



12tt 2 p 2 



e MVfKr ^K(Qi,Qi) (78) 



with 



K(QlQ 2 2 ) = - / J(Qi, z)J{Q 2 , z) d z a{z) dz , (79) 
Jo 

where J(Q,z) is the non-normalizable mode, see 
Eq. ([51]), It satisfies EOM given by Eq. (J30J), and is nor- 
malized by J~(0, z) = 1 for Q — 0. Note that deriving the 
result for K(Q 2 , Q?,) we used the fact that both isoscalar 
and isosinglet vector mesons are described by the same 
EOM and b.c. 

QCD axial anomaly requires: K QCD (0, 0) = 1. In- 
deed, in the present extended holographic model, we get: 



#(0,0) 



d z a(z) dz = a(0) = 1 



(80) 



This result for the anomalous form factor is very similar 
to that obtained in our paper [13] , where we worked out 
a CS extension of the hard- wall model of Refs. @, Q. 
However, in the present model, we do not have a bulk 



K(Q 2 ,0) 



z o Jo 



z J(Q,z)dz . 



(81) 



It is easy to notice that this expression for K(Q 2 ,0) 
coincides with the expression (f6"4")) for F- K (Q 2 ), i.e., the 
anomalous form factor K(Q 2 ,0) in this model coincides 
with the pion EM form factor! 

The slope of the anomalous form factor defined as 



dK{Q 2 ,Q) 



(76) 



dQ 2 



in the present model is given by 



3 2 =0 



16 



^.2 2 

m n z 



(82) 



(83) 



Numerically, we have ~ 0.035 (compare with the re- 
cent result in Ref. [HI). This number is not very far 
from the central values of two last experiments, a T = 
0.026 ± 0.024 ± 0.0048 |39j, a, = 0.025 ± 0.014 ± 0.026 
pol ]. but the experimental errors are rather large. An ear- 
lier experiment [4l[ produced a w = —0.11 ± 0.03 ± 0.08, 
a result whose central value has opposite sign and much 
larger absolute magnitude. In the spacelike region, the 
data are available only for the values Q 2 ^ 0.5 GeV 2 
(CELLO [13) and Q 2 > 1.5 GeV 2 (CLEO |3|) which 
cannot be treated as very small. The CELLO collabo- 
ration [33] gives the value a T = 0.0326 ± 0.0026 that is 
very close to our result. To settle the uncertainty of the 
timclike data (and also on its own grounds), it would 
be interesting to have data on the slope from the space- 
like region of very small Q 2 , which may be obtained by 
modification of the PRIMEX experiment [HJ at JLab. 



C. pcjn coupling 

Substituting V y (x,z) — g^u>^(x)ipi(z) and VHx, z) 
g5P^(x)ipi(z) into Eq. (|T6"j) . we obtain 



-AdS 



Ned 

2tt 2 z 2 U 



dz z ip 1 



Here, we introduced the dimensionful pion field n a = 
firTr a . This lagrangian is similar to that obtained in the 
hidden local symmetries approach [3~ij (see also a review 
[35l]). Thus, we may write that 



N c gj M p 

2lT 2 Zq/tt Jo 

N c g 2 5 M p 

7T 2 U 



d.r 



dz z 
3 J 2 (70,12;) 



(84) 



^i 2 (7o,i) 
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The numerical value of this coupling is g puJ7r ~ —12.77 
(for g$ = V2ir) . The value of this coupling and especially 
its sign have important phenomenological implementa- 
tions, see e.g. Ref. 36]. However, one cannot directly 
measure this coupling constant, since the decay to — » pir 
is energetically forbidden. 

D. Large-Q 2 behavior 



propagators, which is the same (free-field-like) in all the 
models, and by the value of the pion wave function </?(0) 
at the origin, which is fixed by the pion decay constant. 
Numerically, however, the results of the present model 
based on ([85]) would depend on whether one takes the 
experimental value for f n or substitutes gf/ 2 by 2/z 2 . 
Again, the situation is the same as in the pion EM form 
factor case. 



Equation (|HTj) formally gives prediction for the 
K(Q 2 ,0) form factor at all Q 2 , and it is interesting 
to compare it with the monopole fit K CLEO (Q 2 , 0) = 
1/(1 + _Q 2 /A 2 ) (fitted value is = 776 MeV) of CLEO 
data [38| which extend to Q 2 ~ lOGeV 2 . The compari- 
son demonstrating a rather good agreement is shown in 
Fig. 2, where the value zo = 1/(323 MeV) is taken for the 
AdS/QCD curve. 



Q 2 K(Q 2 ,0) (GeV 2 ) 



0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



8 10 



Q 2 (GeV 2 ) 



FIG. 2: Anomalous form factor Q 2 K(Q 2 ,0) from the holo- 
graphic model (solid curve, red online) in comparison with 
the monopole fit Q 2 K Taono (Q 2 , 0) = Q 2 /(l + Q 2 /0.60GeV 2 ) 
of CLEO data (dashed curve, blue online). 

One may also consider the general case of large vir- 
tualities, when Q\ = (1 + uj)Q 2 and Q\ = (1 - ui)Q 2 , 
with to being a fixed parameter, — 1 < u < 1 and Q 2 
large. In this limit, the bulk-to-boundary propagators 
J(Qi,z),J'(Q2 1 z) in Eq. (|79"1) may be substituted by 
their free-field version JC{Qz) = zQKi(Qz). Since JC{Qz) 
exponentially ~ e~® z decreases for large z, the z-integral 
is dominated by small z ~ 1/Q, and Eq. ([79"]) converts 
into 



K{{1 + lu)Q 2 ,{1-lu)Q 2 ) 
2 



(85) 



4Q 2 

Using Eq. (j43|) to substitute the overall factor 2/z$ by 
g 2 f 2 , we obtain exactly the result derived in our earlier 
paper J26JI within the extension of the hard- wall model of 
Refs. [a, Q. This coincidence in analytic form is analo- 
gous to that observed in the pion electromagnetic form 
factor case. Indeed, the large-Q 2 asymptotics is deter- 
mined by the large-Q 2 behavior of the bulk-to-boundary 



V 7 !-^ 2 / d X X 3 K! (xVT+U)K ! ( X VT^J) 



V. SUMMARY AND CONCLUSIONS 

Working within the framework of the holographic dual 
model of QCD proposed by Hirn and Sanz, we study form 
factors of the pion, namely, the electromagnetic form 
factor of charged pions and the anomalous form factor 
7*7*7T° of the neutral pion. In order to calculate the lat- 
ter, we extend the Hirn-Sanz model by incorporating the 
Chern-Simons term into the original 5D action. 

Due to a simple form of the pion "wave function" a(z), 
the pion form factor can be written in explicit analytic 
form involving the modified Bessel function. This ana- 
lytic expression gives a simple formula for the pion charge 
radius in terms of the hard- wall scale zq, which is the 
only parameter of the model fixed by the value of the 
p-meson mass. Written as a function of zq and Q , the 
prediction of the present model is in a good agreement 
with experiment both for low and high Q 2 values. We 
also established that the low energy coupling constant 
9piriT in the present model is in better agreement with 
experiment than the result of the hard- wall AdS/QCD 
model of Ref. @]. 

We extended the model of Hirn and Sanz by adding 
the Chern-Simons term and demonstrated that such an 
extension correctly reproduces QCD anomaly. There 
was no need to introduce an IR counterterm which was 
required in our previous paper [271 ] . where we worked 
in the framework of the hard- wall AdS/QCD model of 
Refs. [6j, la] . We also observed that the anomalous pion 
form factor with one real and one virtual photon with 
momentum transfer Q 2 in the present model is given by 
exactly the same analytic expression as the form factor 
of the charged pion evaluated for the same Q 2 . This out- 
come may be partially due to a very simple form of the 
pion wave function. 

We calculated the Q 2 - s l°P e of the anomalous form fac- 
tor predicted by the pres ent model, which was considered 
earlier in our paper [271 ] and was also recently discussed 
in the Ref. [32j |. In addition, we calculated the value of 
the g P ujTr coupling, which is important for phenomenolog- 
ical considerations. Finally, we showed that in the large 
Q 2 -region we reproduce the same results as in case of the 
hard- wall model of Refs. @, [|[ . 

It is encouraging to establish that such a simple model 
containing just one free parameter, the confinement scale 
Zq, produces the results which are in good agreement with 
experimental findings. Moreover, most of the important 
expressions can be represented analytically without mak- 
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ing any approximations. One can think that the role of 
the scalar field in the hard- wall AdS/QCD model is now 
played by the appropriate b.c. on the IR. However, with 
this simplicity, we loose information about the chiral con- 
densate and the dependence of the observables on it. This 
dependence was studied in our earlier papers [H, H3] • 
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